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Abstract 

We present a search for the decay B~ — > t~V in a 140 fb _1 data sample collected at the 
T(45) resonance with the Belle detector at the KEKB asymmetric B factory. Combinatorial and 
continuum backgrounds are suppressed by selecting a sample of events with one fully reconstructed 
B. The decay products of the other side B in the event are analyzed to search for a B~ — > t~V 
decay. We find no significant evidence for a signal and set a 90% confidence level upper limit of 
Br(B~ — > t~T>) < 2.9 x 10~ 4 . All results are preliminary. 

PACS numbers: 
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The purely leptonic decay B~ — > £~T7 is of particular interest since it provides direct 
measurement of the product of the Cabibbo-Kobayashi-Maskawa(CKM) matrix element V u b 
and the B meson form factor fg. In the Standard Model(SM), the branching fraction of the 
decay B~ — > £~V is given as 

Br{B~ - fr) = f 1 - 4) 2 flW^r B (1) 

87T y m B J 

where Gp is the Fermi coupling constant, mi and m# are the charged lepton and B meson 
masses, tb is the B~ lifetime. The dependence of the lepton mass arises from helicity 
conservation and helicity suppresses the muon and electron channels. 

In the extension of the Standard Model, one expects significant modification to the B~ — > 
t~V decay branching fraction. In the two-Higgs doublet model, the decay can occur via a 
charged Higgs particle. The B~ — > t~V branching fraction is given as 

Br[B- - r-V) = ^± t , _ ^) }l\V^r B x r H , (2) 



r H =[l- , , ,1^ 1 0) 



where r# is defined as 

tan 2 (3 mj 
1 + e~ tan (3 m 2 H t 

and tan/5 is the ratio of vacuum expectation values of two Higgs doublets 0. Once we get 
an upper limit on Br(B~ — > T~v), we can give a constraint on tan/3 and ran. Similarly, 
in i?-parity violating extensions of the MSSM, B~~ — > r~V may be mediated by scalar 
supersymmetric particles^]. Hence, upper limits on the B~~ — > t~V branching fraction 
constrain i?-parity violating couplings. 

No evidence for an enhancement relative to the Standard Model prediction was ob- 
served inprevious experimental studies by CLEO 0,3, ALEPH0, L30, DELPHI and 
BABAR|8J]. The most stringent upper limit has been achieved by the BABAR Collaboration 
: Br{B~ — > t~u) < 4.2 x 10 -4 at 90% C.L. from a sample of fully reconstructed B and 
semi-leptonic decays. 

We use a 140 ftT 1 data sample containing 152.0 x 10 6 B meson pairs collected with the 
Belle detector at the KEKB asymmetric energy e+e" (3.5 on 8 GeV) collider operating 
at the T(4S") resonance (y/s = 10.58 GeV). The Belle detector is a large-solid-angle mag- 
netic spectrometer consisting of a three-layer silicon vertex detector, a 50-layer central drift 
chamber (CDC), a system of aerogel threshold Cerenkov counters (ACC), time-of-flight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals 
(ECL) located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. 
An iron flux-return located outside of the coil is instrumented to identify K\ and muons. 
The detector is described in detail elsewhere III 



Fully reconstructed B mesons, B rec , are observed in the following decay modes: B + — > 
D^°tt + , p + , D^*°af and . D° candidates are reconstructed as D° — > K + ir~, 

K + iT~7r , A' + 7r~7r + 7r~ , K®7i°, K®7t~tt + , K®7t~7t + 7t and K~K + . D*° mesons are recon- 
structed by combining the D° candidates with a pion or a photon. D$ candidates are recon- 
structed in the decay modes Dg — > K S K + and K + K~ii + , and D* s + mesons are reconstructed 
by combining the Dg candidates with a photon. All the tracks and photon candidates in the 
event not used to reconstruct the B rec are studied to search for B~~ — > r~V. The advantage 
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FIG. 1: Distribution of the beam energy constrained mass in data for fully reconstructed B mesons 
(histogram). The solid curve shows the result of the fit and the background components as the 
dotted curve. 

of having a sample of fully reconstructed B meson is to provide a strong suppression of the 
combinatorial and continuum background events. The disadvantage is the low efficiency of 
full B meson reconstruction (about 0.3%). Charged B pair events are generated from T(4s) 
resonance (y/s ~ 10.58 GeV), where the B + or B~ is generated with specific momentum and 
energy. Selection of the fully reconstructed B candidates is made according to the values 
of two variables : the beam constraint mass M bc = J E% eam — p B and the energy difference 
AE = E B — E beam . Here, E B and pb are the reconstructed energy and momentum of the 
fully reconstructed B candidate in the center-of-mass (CM) system, and E beam is the beam 
energy in the CM frame. 

The Mhr distribution of reconstructed B candidates is fit with the sum of an Argus 
function ^lj and a Crystal Ball function 0). The Argus function models the continuum and 
combinatorial background whereas the Crystal Ball function models the signal component, 
which peaks at the B mass. The purity is denned as S/(S + B), where S (B) is the 
number of signal (background) events for M bc > 5.27 GeV/c 2 , as determined from a fit. 
Figure ^ shows the M bc distribution for all B rec candidates in our data set. The yield 
Nb+ b- of the sample containing one B rec is determined as the area of the fitted Crystal 
Ball function. We obtain N b +b- = (2.40 ± 0.15) x 10 5 and 0.57 of the purity, where the 
uncertainty on N b +b- is dominated by systematic errors. We define the B rec signal region 
to be —0.08 < AE < 0.06 GeV and M bc > 5.27 GeV/c 2 . In order to avoid experimenter 
bias, the signal region in data is blinded until the selection criteria are finalized. 

Once a reconstructed B candidate has been identified, B~ — > t~V signal candidates are 
selected by considering all tracks and clusters in the event which are not used in the fully 
reconstructed B. Candidate events are required to have one or three signal-side charged 
track(s) with the total charge which is opposite that of the reconstructed B. In the events 
where a B rec is reconstructed, we search for decays into a r plus a neutrino. The r lepton 
is identified in the following decay channels: t~ — > fi~uu, r~ — > e~vi>, t~ — > tt^u, t~ — > 
7r _ 7r°i/, and r~ — > ti~tt + ti~v. 

We require the charged particles to be identified as leptons or pions. For the lepton and 
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single-pion modes we reject events with ir° mesons in the recoil against B rec . The event 
is required to have zero charge and Eecl less than a certain value (Eecl < 1.0 GeV for 
t~ — > \i~vv and e~vv y Eecl < 1.2 GeV for t~ — > n~v and 7r~7r + 7r~z/ and Eecl < 2.2 GeV 
for r" -> 7r"7r°z/) where E SCL is defined as E SCL = E tot - E Brec - E track , the energy 
deposition in the ECL calorimeter. An additional requirement E n o/E E cL > 0.2 is applied 
for t~ — > 7T~7r°z/. Further requirements are made on the total momentum of the track(s) in 
the CM frame (p w - > 1.0 GeV/c for r~ — > n^u, p 7r - 7r + 7r - > 1.2 GeV/c for r~ — > ix~tx+ix~v), 
the total missing momentum of the event {P > 1.0 GeV/c for r~ — > 7r~z/ and 7r~7r ^, 
^ > 1.2 GeV/c for t~ — > 7r~7r + 7r _ z/) and the invariant mass of two or three pions 0.55 < 
< 0.95 GeV/c 2 and 1.0 < < 1.4 GeV/c 2 . The selection efficiencies for the r 

decay channels we consider in this analysis are determined from Monte Carlo simulations of 
B~ — > t~v events. We compute the efficiency as the ratio of the number of events surviving 
each of our selections over the number of fully reconstructed B. 

The expected background is composed of events from continuum and combinatorial back- 
ground, along with fully reconstructed B meson events. Backgrounds consists primarily of 
B + B~ events in which the fully reconstructed B has been correctly reconstructed. This 
"peaking" background is directly determined from Monte Carlo simulations of B + B~ events. 
The continuum and combinatorial background is determined from the number of events in 
Mb c sideband data, scaled by the ratio of the areas of the fitted Argus function in the signal 
and sideband regions. We fit Mf, c distributions after preselection and assume the ratio of 
the fitted Argus is unchanged after all selection criteria have been applied. This assumption 
is consistent with the observed distributions in the Monte Carlo as well as the data in AE 
sideband. 

The main sources of uncertainty we consider in the determination of the Br(B~ — > 
t~V) are uncertainty in the number of B + B~ events with one reconstructed B, uncertainty 
in the determination of the signal efficiency and uncertainty in the determination of the 
number of expected background events. The number of B + B~ events is determined as 
the area of the Crystal Ball function fitted to the Mf, c distribution. Using a Gaussian 
function as an alternative fitting function, we obtain a relative change in the number of 
events and this difference is assigned as the systematic uncertainty on the number of B + B~ 
events. A systematic error due to uncertainty in the amount of neutral B background in the 
number of B + B~ is also considered. The main contribution to the systematic uncertainties 
in the determination of the efficiencies come from uncertainty on tracking efficiency, Monte 
Carlo statistics, Eecl energy and particle identification. The uncertainty in the expected 
background comes from Monte Carlo statistics, Eecl energy uncertainty and use of the 
Argus fit function. 

Figure 121 shows the Eecl distributions in the data after all selection requirements ex- 
cept the one on E E cl have been applied compared with the expected background. Each 
distribution refers to a different selections and the plots show no evidence of signal in data. 
We find a total of 28 candidates in the signal region where 33.5 ± 5.0 background events 
are expected. This uncertainty in the background expectation is due to both statistical and 
systematic errors. 

In order to extract the upper limit on the branching fraction for B~ — > t~V, we combine 
the results of the different selections. We use the likelihood ratio estimator, Q. Here, Q is 
defined as C(s + b)/C(b), where C(s + b) and C{b) are the likelihood functions for signal plus 
background and background only hypotheses, respectively. The likelihood functions C(s + b) 
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FIG. 2: Eecl distributions in the data after all selection requirements except the one on Eecl 
have been applied. The vertical arrow is the requirement on the Eecl m each selection. 



and C(b) are denned as 



t=l ' b i- 



n ch —biUni 



i=l 



(4) 



where n c h is the number of selection channels, Sj and hi are the expected number of signal and 
background events, respectively, and rij is the number of observed events in each channel. 
The number of signal events Sj can be written ■ N B + B - ■ Br(B — > r i/), where 

£j is the selection efficiency for the i-th channel and Nb+b- is the number of B + B— events 
with one reconstructed B rec . We set a 90% CL. upper limit using a simple Monte Carlo 
generating random experiments for different values of the branching fraction Br(B~ — > t~v). 
The confidence level for the signal hypothesis can be written as 



CL 



s+b 



Nr 



CL h 



Qb<Q 



(5) 



where NQ a+b<Q and NQ b<Q are the number of the generated experiments which have a likeli- 
hood ratio less than or equal to the measured one, in the signal plus background hypothesis 
and background only hypothesis, respectively. The 90% CL. upper limit is obtained for 
CL S = 1 — 0.9, where Gaussian systematic uncertainties are incorporated. This "Modified 
Frequentist approach" or "CLs method" method is described in detail in reference [l3j |. 
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FIG. 3: The confidence level for the signal hypothesis CL S as a function of Br(B~ — > t~v) 
obtained from the likelihood approach. The solid curves are the observed results, the dashed curves 
the expected. The shaded areas represent the symmetric la and 2a bands. The intersections of 
the curves with the horizontal line at CL S = 0.1 give the limits on Br(B~ — ► t~v) at the 90% 
confidence level. 



Figure El shows the confidence level for the signal hypothesis CL S as a function of 
Br(B~ — > t~v) obtained from the likelihood approach. The solid curves are the observed 
results, the dashed curves the expected. The shaded areas represent the symmetric la and 
2cr bands. Intersection of the curves with the horizontal line at CL S = 0.1 gives the limits 
on Br(B~ — > t~v) at the 90% confidence level, with systematic uncertainties taken into 
account. The limits we observe lie within the la expected regions. 

Table U] shows the branching fraction times efficiency for signal, expected background, 
observed events and upper limit on the branching fraction for each r decay mode. We 
obtain a combined limit on the branching fraction of 

Br(B~ -> t~V) < 2.9 x 10 -4 at the 90% C.L. (6) 

In the two-Higgs doublet model, the branching fraction Br(B~ — > t~v) is enhanced by a 
factor of [1 — (m B /m H ) 2 tan 2 /3} 2 for e = in equation Q. With m B = 5279 MeV/c 2 and 
Br(B~ — > t~V) = 7.5 x 10~ 4 from the Standard Model prediction, we set the constraint 

tan p nnn /n ir , 2 



M H ± 



< 0.33 (GeV/c 2 )- 1 (7) 



from the obtained limits on the branching fraction of B~ — > t~V. Figure |U shows the 90% 
C.L. exclusion boundaries in the \M H + , tan 0} plane obtained from J7J) compared with other 
experimental searches at LEP W| and at the Tevatron(l5|. The direct search for charged 
Higgs bosons at LEP gives the constraint M#± > 78.6 GeV/c 2 . Figure El shows the summary 
of searches for B~ — * t~V and upper limits obtained by other experiments compared with 
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t Decay Mode 


Efficiency xBR(%) 


Background Observed 
Expected Events 


UU&cl VcU lllilll 

^yu /o\-j.i-i. j 


T~ — > H~VU 


9.2 ± 1.5 


9.8 ±2.9 


6 


3.0 x 10~ 4 


t~ — > e~vv 


O.O ±1.0 


9.4 ±2.9 


10 


a a v> 1 n — 4 
4.0 X iU 


T~ — > -K~VV 


4.1 ±0.4 


5.4 ±2.1 


6 


7 9x1 n -4 


T~ — > 7l~7V°V 


1.8 ±0.2 


4.1 ± 1.6 


3 


10.5 X 10~ 4 


T~ — > 7T + 7r~7r + Z/ 


1.6 ±0.2 


4.8 ± 1.6 


3 


11.7 x 10" 4 


Combined 


2.9 x 10" 4 



TABLE I: Branching fraction times efficiency for signal, expected background, observed events and 
upper limit on the branching fraction for each r decay mode. The uncertainties on the efficiencies 
and the expected background are due to the combination of limited statistics and systematic errors. 



the corresponding SM prediction and the branching fraction predicted with the charged 
Higgs boson in two parameter sets. Our result is the most stringent on this process. 

In conclusion, we have performed a search for the B~ — > t~V decay in a fully reconstructed 
B sample. The analysis uses the following r decay channels: t~ — > fi~i>v, t~ — > e~vv, 
t~ — > 7r~b>, r~ — > 7r~7r°z/, and t~ — > tt~7t + tt~v. The results of the search in the different 
channels have been combined. No signal is observed and an upper limit has been set : 

Br(B- -> t~V) < 2.9 x 10~ 4 (90% C.L.), (8) 

which represents the most stringent upper limit on this process to date. 
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FIG. 4: The 90% C.L. exclusion boundaries in the [M#+,tan/3] plane obtained from the observed 
upper limit on Br(B~ — ► t~v). 
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